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Abstract 

Protein isolates of pumpl<in {Cucurbita pepo L) seeds were hydrolyzed by acid protease to prepare antioxidative peptides. 
The hydrolysis conditions were optimized through Box-Behnken experimental design combined with response surface 
method (RSM). The second-order model, developed for the DPPH radical scavenging activity of pumpkin seed hydrolysates, 
showed good fit with the experiment data with a high value of coefficient of determination (0.991 8). The optimal hydrolysis 
conditions were determined as follows: hydrolyzing temperature 50°C, pH 2.5, enzyme amount 6000 U/g, substrate 
concentration 0.05 g/ml and hydrolyzing time 5 h. Under the above conditions, the scavenging activity of DPPH radical was 
as high as 92.82%. 
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Introduction 

With the improvement of people's living standard, many kinds 
of diseases caused by reactive oxygen species (ROS) and free 
radicals are bringing serious damage to human health [1,2]. In 
recent years, preparation of antioxidative compounds has become 
a research hotspot in biology, medicine and food science [3,4]. 
Nowadays, more and more attention has been focused on plant 
protein hydrolysates with antioxidative properties [5,6], which are 
considered as potential natural substitutes of artificial antioxidants 
(BHA, BHT, TBHQ), as well as exogenous antioxidants in human 
nutrition [7]. 

Considering the nutritional quality and cost, many legumes, 
oilseeds and their by-products of the oil industry are suitable as 
protein sources to produce protein hydrolysates for food or non- 
food application. Some materials have been extensively studied 
including soybean protein [8], peanut protein [9,10], cottonseed 
protein [1 1], etc. Many studies show that the antioxidative activity 
of individual peptide depends on amino acids in sequence and 
molecular mass [12-14]. However, the hydrolysate with antiox- 
idative properties of one kind of protein may also be affected by 
some other factors, such as pH, temperature, the enzyme amount 
and hydrolysis time. When many different factors need to be 
considered, the response surface method (RSM) is an effective tool 
for optimization of conditions [15,16]. 

Pumpkin seed is one of the most important oilseeds in many 
countries (India, China, etc.), which is also a popular snack. Nearly 
90% of the total output in these countries is used for extraction of 
oil [17], leaving a large amount of residue which contains 60%- 
65% of high value protein. Although some defatted flour is utilized 
increasingly, a huge amount of the residue is still being discarded 
causing waste and pollution. Additionally, pumpkin seeds could 
also be used as raw material for production of high quality protein 



products which are nutrition supplements for food formulations, as 
well as functional agents [18-21]. 

Most notably, pumpkin seed is a valuable raw material for 
preparation of antioxidative peptides owing to its high protein 
content, availability and low cost. In this paper, antioxidative 
peptides were prepared by hydrolysis from pumpkin seed. The 
effects of enzyme amount, substrate concentration and hydrolysis 
time on the DPPH radical scavenging ability of pumpkin seed 
hydrolysates were investigated. Further more, the optimal 
hydrolysis conditions were established through RSM. 

Materials and Methods 

Materials 

Pumpkin seeds [Cucurbita pepo L.) were provided by "JinlUai" 
Co., Shanxi, China. Enzyme was acid protease with a declared 
activity of 100,000 U/g purchased from Imperial Jade Bio- 
technology Co., China. l,l-diphenyl-2-picrylhydrazyl (DPPH) 
was obtained from Sigma Chemicals Co. AH the other chemicals 
are chemically pure (chemical reagent). 

Preparation of Protein Isolate 

The protein isolate was prepared in alkaline solution followed 
by isoelectric precipitation, as described by Sathe [22] with some 
modifications. Pumpkin seeds were dehuUed to get kernels which 
were ground in a hammer mUl. The ground kernels were then 
defatted by supercritical carbon dioxide to give the defatted flour, 
which was allowed to mix with water in the ratio of 1/30 (w/v). 
The pH of the mixed solution was adjusted to pH 11 by IN 
sodium hydroxide (NaOH). The mixed solution was stirred at 
50°C for 1.5 h, centrifuged at 5000 r/min for 20 min. The 
supernatant was collected and adjusted to pH 5.3 by IN 
hydrochloric acid (HCl). Then the generated pH 5.3 solution 
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Table 1. Factors and levels in the response surface design. 



Factor Symbol Levels 

-1 0 1 

Enzyme amount (U/g) x, 2000 5000 8000 

Substrate concentration (g/mi) 0.02 0.04 0.06 

Time (h) Xj 2 4 6 

doi:l 0.1 371 /journal.pone.0092335.t001 

was centrifuged at 5000 r/min for 20 miii. The supernatant was 
discarded and the precipitate was collected and dried using a 
freeze drier. 

Experimental Design 

In preliminary studies, amount of enzyme, substrate concen- 
tration and time were found to have a significant influence on the 
degree of hydrolysis (DH) and DPPH radical-scavenging activity. 
The effect of three independent variables on the scavenging 
activity and their optimum levels were investigated by RSM, using 
Box-Behnken experimental design. Here, 15 experiments were 
carried out (with three replicates at the centre of the design). The 
range and central point values of the independent variables 
(Table 1) were assorted based on the results of preliminary 
experiments. The full experimental design with respect to 
predicated values and experimental values for the response (DPPH 
radical-scavenging activity) is presented in Table 2. 

Hydrolysis of Protein Isolates 

The enzymatic hydrolysis was carried out as stated in Table 1, 
according to the experimental plan in Table 2. The hydrolysis was 
performed in a glass vessel under controlled conditions. The 
mixture was incubated at 50°C and pH 2.5. Then the mixture was 



Table 2. Box-Behnken design matrix and the response values for the DPPH free radical scavenging ability (%) of pumpkin seed 
hydrolysates. 



Treat 


Variable levels 




DPPH free radical scavenging ability (%) 




Xl 






Experimental (Yo) 


Predicted (Yi) 


1 


0 


-1 


-1 


80.05 


79.8937 


2 


0 


-1 


1 


84.47 


85.0288 


3 


0 


0 


0 


91.32 


90.9767 


4 


0 


0 


0 


90.95 


90.9767 


5 


0 


1 


-1 


90.84 


90.2813 


6 


0 


0 


0 


90.66 


90.9767 


7 


-1 


-1 


0 


81.63 


81.6350 


8 


1 


-1 


0 


84.72 


84.3125 


9 


1 


0 


1 


91.51 


91.3588 


10 


-1 


1 


0 


90.78 


91.1875 


11 


1 


1 


0 


92.56 


92.5550 


12 


-1 


0 


1 


90.74 


90.1763 


13 


0 


1 


1 


92.28 


92.4363 


14 


1 


0 


-1 


87.99 


88.5538 


15 


-1 


0 


-1 


85.54 


85.6913 
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separated to 3 portions, which are allowed to heat up to 100°C in 
5 min and react for 2 h, 4 h and 6 h, respectively. 

The mixture with inactivated enzyme was centrifuged at 
5000 r/min for 10 min and the supernatant was used for further 
analysis. 

Determination of the Degree of Hydrolysis 

The degree of hydrolysis was determined according to the 
method of Chalamaiah, M. et al. [23]. An aliquot (10 ml) was 
taken out from the supernatant of protein hydrolysis solution, then 
mixed with 10 ml of 20% trichloroacetic acid (TCA) and 
centrifuged at 5000 rpm for 30 min at room temperature. The 
supernatant was decanted and nitrogen analysis was carried out 
with the micro-Kjeldahl method [24] . The degree of hydrolysis 
(%) can be calculated according to equation (1). 



^^TT«/^ 10%TCA - solublenitrogeninsample 

\Jn\/o)= — -, : X iUU lil 

totalmtrogenmsample 



Scavenging Activity of DPPH Radical 

The scavenging activity of DPPH radical was measured 
according to the method of Shimada [25] with some modifica- 
tions. In each test, an aliquot (2 ml) of sample solution with 
identical concentration was allowed to mix with 2 ml of DPPH 
(0.2 mmol/L in ethanol). The reaction mixture was incubated for 
50 min in dark at room temperature. Absorbance of the resulting 
solution at 5 1 7 nm was measured by UV-Vis spectrometer (UV- 
2800AH, Shanghai, China). Ethanol was used to controlling 
solvent. The radical-scavenging capacity of samples was evaluated 
by increasing/ decreasing the concentration of DPPH radical 
calculated by equation (2). 
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Table 3. Significant test for each regression coefficient of the fitted regression model. 



Variables 


Regression coefficient 


Standard error 


r-value 


P-value 


Significance level 


Constant 


90.9767 


0.3564 


255.244 


<0.0001 


*# 


Xl 


1.0113 


0.2183 


4.633 


0.006 


** 


X2 


4.4488 


0.2183 


20.382 


< 0.0001 


»* 


X3 


1.8225 


0.2183 


8.350 


< 0.0001 


»* 


Xl^ 


-0.7596 


0.3213 


-2.364 


0.064 




x.^ 


-2.7946 


0.3213 


-8.698 


< 0.0001 


** 


X3^ 


-1.2721 


0.3213 


-3.959 


0.011 




Xl X2 


-0.3275 


0.3087 


-1.061 


0.337 




Xl X3 


-0.4200 


0.3087 


-1.361 


0.232 




X2 X3 


-0.7450 


0.3087 


-2.414 


0.061 





•P<0.05, ••P<0.01. 

doi:l 0.1 371/journal.pone.0092335.t003 



Scavengingactivity(%) = 



.'^sample 



X 100 



(2) 



polynomial equation was presented as surface plots in order to 
visualize the relationship between the response and experimental 
levels of each factor and search out the optimum conditions. The 
analysis software used for this study was Minitab 15. (Minitab 
Inc.). 



Statistical Analysis 

Data from the Box-Behnken experimental design (Table 2) was 
analyzed to determine the optimum combination of variables. The 
polynomial model proposed is given as equation (3): 



Y = bo + biXi 



- b2X2 + b3X3 + b4Xi + hsxj 



+ b6X3 +b7XiX2 + bgXiXi + b<)X2X3 



(3) 



where Y is the predicted response variable (scavenging capacity of 
DPPH radical), bo is the regression coefficient of intercept term, 
bi, b2, bs are hnear regression coefficients, b4, bs, be are squared 
regression coefficients and by, bg, bg are interaction regression 
coefficients. 

The significance of each coefficient of the resulted model was 
determined by using the Student t-test and p-value. The 
proportion of the variance expressed by the models is obtained 
by determination of multiple coefficient R^. Finally, the fitted 



Results and Discussion 

Model Fitting 

RSM was apphed to obtain the regression equation, which 
represents an empirical relationship between response (scavenging 
activity of DPPH radical) and the test variables (enzyme amount, 
substrate concentration and hydrolysis time), as given in the 
equation (4). 



Y = 90.9767 +1.011 34xi + 4.4488x2 + 1 .8225x3 - 0.7596xf 

(4) 

-2.7946X?- 1.2721x2 -0.3275x1X2 -0.4200x1X3 -0.7450x2X3 



The significance of each coefficient was determined by using the 
/-values and P-values shown in Table 3. It can be seen that the 
largest effect was caused by the linear term including enzyme 
amount (xj), substrate concentration (X2) and time (x^), followed by 
the quadratic term of substrate concentration (x.^^) and the 
quadratic term of time (x^^). The interaction terms did not have 



Table 4. Analysis of variances for the fitted regression model. 





Source 


df 


Sum of squares 


Mean squares 


F-value 


P-value 


Significance level 


Regression 


9 


230.247 


25.5830 


67.12 


<0.0001 


** 


Linear 


3 


193.084 


64.3614 


168.87 


<0.0001 


** 


Square 


3 


33.808 


11.2694 


29.57 


0.001 


»» 


Interaction 


3 


3.355 


1.1182 


22.9 


0.138 




Lack of fit 


3 


1.687 


0.5623 


5.14 


0.167 




Pure error 


2 


0.219 


0.1094 








Total 


14 


232.152 






R^ = 0.9918 





*P<0.05, **P<0.01. 
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PLOS ONE I www.plosone.org 



3 



March 2014 | Volume 9 | Issue 3 | e92335 



Preparation of Antioxidative Peptides Using RSIVl 




Figure 1. Effect of enzyme amount and substrate concentration. Response surface plots (A) and contour plots (B) of the DPPH radical 
scavenging activity of pumpkin seed hydrolysates affected by enzyme amount and substrate concentration. 
doi:10.1371/journal.pone.0092335.g001 



significant influence (P>0.05), which means that the interaction 
between the different factors did not influence the response. 

The calculated values by using the regression model were 
allowed to do a comparison with experimental values. The 
coefficient of determination of calculated values was 0.9918. 
This indicates that the fitted model could explain 99.18% of the 
total variability within the range of values studied. The analysis of 
variance (Table 4) showed P-value was less than 0.01, which 
implied that the model itself is extremely significant. Meanwhile, 
the statistical analysis data revealed that linear and quadratic terms 
were extremely significant (P<0.01) while the interaction coeffi- 
cients not (P>0.05). The lack of fit was also non-significant (P> 
0.05), which further validated the model. 

Optimization of the Process 

The three-dimensional response surface graphs and contour 
maps were shown in Fig. 1 to Fig. 3, which illustrate the interactive 
effects of the independent variables on the scavenging capacity of 
DPPH radical. The shape of contour reflects the strength of the 
interaction effects, circular means interaction is not significant, 
while oval is opposite. 

Fig. 1 showed the effect of enzyme amount and substrate 
concentration on the scavenging capacity of DPPH radical, we 



learn that substrate concentration had greater influence on the 
response. Higher values of scavenging capacity were observed 
when the enzyme amount was 4000 to 7000 (U/g) and the 
substrate concentration was 0.04 to 0.06 (g/ml). Fig. 2 showed the 
effect of enzyme amount and time. Scavenging capacity of DPPH 
radical boosted with the increase of enzyme amount and time. 
Fig. 3 demonstrated the effect of substrate concentration and time. 
When substrate concentration was 0.04 to 0.06 g/ ml and time was 
3 to 6 h scavenging capacity of DPPH radical exhibited higher 
values of 92% to 96%. 

These results implied that the response surface had a maximum 
point in the experimental range of the independent variables. The 
precise coordinates of optimum (i.e the values of three indepen- 
dent variables) were obtained by analytical procedure: enzyme 
amount 6181.82 U/g, substrate concentration 0.0543 g/ml, time 
4.87 h. Thus, we can learn that the maximum predicted 
scavenging capacity under the optimal condition was 93.17%. 

Confirdmatory Tests 

To confirm the validity of the suggested mathematical model, 
the trial experiments were conducted under the predicted 
optimal conditions. Take convenience into account, the actual 
optimum experimental parameters were modified as follows: 




Figure 2. Effect of enzyme amount and hydrolysis time. Response surface plots (A) and contour plots (B) of the DPPH radical scavenging 
activity of pumpkin seed hydrolysates affected by enzyme amount and hydrolysis time. 
doi:1 0.1 371/journal.pone.0092335.g002 
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Figure 3. Effect of substrate concentration and hydrolysis time. Response surface plots (A) and contour plots (Bj of the DPPH radical 
scavenging activity of pumpkin seed hydrolysates affected by substrate concentration and hydrolysis time. 
doi:1 0.1 371/journal.pone.0092335.g003 



enzyme amount 6000 U/g, substrate concentration 0.05 g/ml, 
time 5 h, with the temperature of 50°C and pH 2.5. Three 
parallel experiments were performed, and the values of DPPH 
radical scavenging activity were 92.15%, 93.08% and 93.23%, 
respectively. The mean of three replicate determinations was 92. 
82%. The experimental value was in high agreement with the 
predicted value, which showed the validity of this response 
model. 

Conclusions and Discussion 

To conclude, we carried out hydrolysis experiment for 
pumpkin seed protein isolate under 15 conditions by different 
combinations of enzyme amount, substrate concentration and 
hydrolysis time (independent variables). The second-order model 
developed for the DPPH radical scavenging activity of pumpkin 



seed hydrolysates was extremely significant (P<0.01) with a high 
value of coefficient of determination (0.9918). The surface and 
contour graphs indicated that maximum scavenging capacity of 
DPPH radical was 92.82%, which obtained under optimum 
hydrolysis condition (enzyme amount 6000 U/g, substrate 
concentration 0.05 g/ml, hydrolysis time 5 h, t temperature 
50°C and pH 2.5), which could be validated by experiment. 
Thus, we believe that the optimum hydrolysis conditions can be 
obtained by RSM precisely. 
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